Equations of motion approach to decoherence and current noise in ballistic 
interferometers coupled to a quantum bath 
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We present a technique for treating many particles moving inside a ballistic interferometer, under 
the influence of a quantum- mechanical environment (phonons, photons, Nyquist noise etc.). Our 
approach is based on solving the coupled Heisenberg equations of motion of the many-particle 
system and the bath and is inspired by the quantum Langevin method known for the Caldeira 
Leggett model. It allows to study decoherence and the influence of the bath on other properties 
of the interferometer. As a first application, we treat a fermionic Mach-Zehnder interferometer. 
In particular, we discuss the dephasing rate and present full analytical expressions for the leading 
corrections to the current noise, brought about by the coupling to the quantum bath. 

PACS numbers: 73.23.-b, 72.70.+m, 03.65.Yz 



I. INTRODUCTION 

Decoherence, the destruction of quantum-mechanical 
phase coherence by a fluctuating environment, plays an 
important role, ranging from fundamental questions such 
as the quantum-classical correspondence to potential ap- 
plications (like quantum information). For most of the 
recent two decades, the focus of research has been on 
quantum-dissipative systems with few degrees of free- 
dom: the most prominent examples are the single particle 
(e.g. Caldeira-Leggett modeU*^) or a single two-level sys- 
tem (spin-boson model) and other impurity models (e.g. 
the Kondo model in transport through quantum dots). 
However, such a description is no longer adequate when it 
comes to transport interference effects both in disordered 
systems (weak localization, universal conductance fluc- 
tuations) or man-made interference devices (Aharonov- 
Bohm rings, double quantum dot interferometers, atom 
chip interference setups etc.). In those cases, we are deal- 
ing with a many-particle system. As long as this is cou- 
pled to classical noise, we can still use the single particle 
picture. Both the technical effort and the pliysical ideas 
expand considerably when going over to a full quantum 
bath. Up to now, there have been comparatively few 
treatments of quantum-dissipative many-particle systems 
(for examples see^i^i^iSii and references therein). 

In this article, we describe an equations of motion ap- 
proach for ballistic interferometers coupled to a quantum 
bath (Fig. QJ. It is physically transparent, more effi- 
cient than generic methods (like Keldysh diagrams) and 
straightforwardly keeps important physics, such as the 
effects of Pauli blocking in fermionic systems. It may 
be applied to describe decoherence (or dephasing; we use 
the terms interchangeably) and, in general, to calculate 
the current noise and other higher-order correlators of 
the particle field. 

We have already introduced this method in a re- 
cent short article^, and applied it to the electronic 
Mach-Zehnder interferometer realized at the Weizmann 
institute^ , discussing the loss of visibility in the current 
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Figure 1: (a) The Caldeira-Leggett model (single particle and 
oscillator bath) and (b) a ballistic many-particle system sub- 
ject to a quantum noise potential V{x,t). 



interference pattern and (briefly) the effects on the cur- 
rent noise. The purpose of the present paper is four- fold: 
(i) to relate our many-particle method to the Quantum 
Langevin equation as it is known for a single particle in 
the context of the Caldeira-Leggett model (Sec. Ill A|l . 
(ii) to provide more details of the method and on how 
to evaluate the resulting expressions using perturbation 
theory (Sec. IIV C|) . (iii) to derive and present full analyt- 
ical expressions for the leading correction to the current 
noise of a MZ setup coupled to a quantum bath (Sec. 
and (iv) to add to our previous brief discussion of 
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the current noise (Sec. IIV E|l . 



II. EQUATIONS OF MOTION APPROACH TO 

DECOHERENCE IN BALLISTIC 

INTERFEROMETERS 

A. Brief reminder of the quantum Langevin 
equation 

The quantum Langevin equation can be employed to 
solve the Caldeira-Leggett modeUiS of a single particle 
coupled to a bath of harmonic oscillators. Briefly, the 
idea is the following, when formulated on the level of 
Heisenberg equations (where it is formally exact). The 
total quantum force F acting on the given particle, due 
to the bath particles, can be decomposed into two parts: 



i^(t)=F(o)(t)+ / D^{t^t')x{t')dt' (1) 

J —OQ 

The first describes tlie intrinsic fluctuations, present even 
in absence of the couphng. It derives from the solution to 
the free equations of motion of the bath oscillators, with 
fluctuations due to the stochastic initial conditions. For 
example, the force might be a linear superposition of nor- 
mal oscillator coordinates, i^(o)(0 — J2j 9jQj{o){t) with 
Qj(0)W - Qmi^)cos{n,t) + (P,(o)(0)/M,l],)sin(f^,i). 
The fluctuations of Qj(q){Q) and Pj(-o)(0) includes both 
thermal and quantum (zero-point) fluctuations. The sec- 
ond part of the force is due to the response of the bath to 
the particle's motion (here: the i-coordinate, if the cou- 
pling is of the type Fx). We will call it the "back-action" 
term, and it gives rise to features such as mass renormal- 
ization and friction. As the bath equations of motion 
are linear (since we are dealing with a harmonic oscil- 
lator bath), the response is linear for arbitrary coupling 
strength, and the resulting equation ^ for the force is 
valid on the operator level (not only for averages) . In this 
way, one has "integrated out" the bath by solving for its 
motion. Plugging the force F into the right-hand-side 
(rhs) of the Heisenberg equation of motion for x yields 
the quantum Langevin equation: 



linearized dispersion relation, as this is sufficient to de- 
scribe decoherence (neglecting acceleration/retardation 
effects). We start from Heisenberg's equations of motion 
for the particles and the bath. A particle field moving 
baUistically at constant speed (see Fig^ (b)) obeys the 
following equation: 



i{dt — VFdx)iip{x, t) = I dx' K{x — x')V{x\ t)^{x' , t) 



(4) 

where V evolves in presence of the interaction, see be- 
low. Here vp would be the Fermi velocity in the case 
of fermions, or the velocity with which bosons have 
been injected into the interferometer (e.g. as a BEC 
cloud in an atom chip, or the speed of light for pho- 
tons). We must consider states within a finite band, 
thus K{x — x') ~ {'ip{x)^'ip'^ {x')} / 5{x — x') (written 
for fermions, analogous for bosons). Nevertheless, for 
the purpose of our subsequent leading-order approxima- 
tion, it turns out we can replace the right-hand side by 
V{x,t)%l)[x,t) (neglecting, e.g., velocity-renormalization 
in higher orders) . The corresponding formal solution de- 
scribes the accumulation of a random "quantum phase" : 



d'^xjt) 
' dt^ 



= F{t)-U'{x{t)) 



(2) 



In practice, this equation can only be solved for a har- 
monic potential U{x), i.e. for a free particle or a har- 
monic oscillator. This is why the range of applications 
of the quantum Langevin equation is usually rather re- 
stricted. For the example of a harmonic oscillator (bare 
frequency wo), we have, with the help of Eq. |^ and 
after going to frequency space: 



?/'(a;,i) — Texp 

ip{x -VF{t-to),to) 



i j dtiV{x ~~ vpit — ti),ti) 

to 



(5) 



In contrast to the case of classical noiseiS, the field V 
contains the response to the particle density, in addition 
to the homogeneous solution V(o) of the equations of mo- 
tion (i.e. the free fluctuations): 



x{lli) 



F{0}{^) 



m{uj^ - uj^) - D^icu) 



(3) 



One can obtain averages of moments of x (and p) by 
plugging in the formal solution and employing the corre- 
lator of F(o) (t) (using the Wick theorem for higher-order 
correlators). 

In the case of a many-particle system, it is the density 
h{x) = tp^x)ip{x) that couples to a scalar noise potential 
V{x). The place of x and F in the quantum Langevin 
equation for a single particle is thus taken by the particle 
field ■0 and V, respectively. 



B. Coupled equations for the many-particle system 
and the bath 



Let us now turn to the case of many particles (fermions 
or bosons) traveling baUistically inside the arm of an in- 
terferometer. We will assume chiral motion and use a 



Vix,t)^V^o)ix,t)- 



dt' D'^{x,t,x',t')n{x',t'). (6) 



Here D^ is the unperturbed retarded bath Green's func- 
tion, £'^(1,2) = -ie{ti -t2)([V{l),Vi2)]y where V- 

correlators refer to the free field. This (exact) step 
is analogous to the derivation of an operator quantum 
Langevin equation, see above. Together with ([SJ, it cor- 
rectly reproduces results from lowest-order diagrammatic 
perturbation theory. 

Below we will apply our approach to the fermionic 
Mach-Zehnder interferometer, though the influence of 
quantum noise on bosonic interferometers (like inii) rep- 
resents another interesting future avenue of research. We 
note that recently a different kind of quantum Langevin 
method has been developed for transport through quan- 
tum dotsi^. 



0^ 




Figure 2: Schematic of the Mach-Zehnder setup, with beam 
sphtters A,B, input ports 1,2, and output ports 3,4. 



III. APPLICATION TO THE MACH-ZEHNDER 
INTERFEROMETER 

A. Introduction 



observed features such as the increase of the dephasing 
rate with rising bias voltage could not be studied, as this 
is a true many-body effect (see below). 



B. The model 

We consider a model of spinpolarized fermions, mov- 
ing chirally and without backscattering through an in- 
terferometer at constant speed vp (see Fig. 121. The two 
beamsplitters A and B connect the fermion fields 4'a of 
the input (a — 1,2) and output {a = 3,4) channels to 
those of the left and right arm {a — L,R), which we take 
to be of equal length I: 



Mo,t) = rAMO,t) + tAMO,t) (7) 

^i?(0,i) = tAMO,t)+rAMO,t) (8) 

Mht) = rBe''^i^L{l,t) + tBi^Ril,t) (9) 

Ml,t) = tBe"^Ml,t) + rBMl,t) (10) 

The transmission (reflection) amplitudes t^/B (^a/b) f^l" 
fill t*rj = —tjT* due to unitarity, and we have included 
the Aharonov-Bohm phase difference (j). The input fields 
a — \,2 are described by Fermi distributions /i_2, where 
the chemical potential difference defines the transport 
voltage: eV = ^x — ^j.2- We have 



In contrast to the usual mesoscopic Aharonov-Bohm 
ring setups, the recently realized Mach-Zchnder interfer- 
ometer for electronsSii^ offers the possibility to study an 
ideal two-way interference geometry, with chiral single- 
channel transport and in the absence of backscattering. 
The loss of visibility with increasing bias voltage or tem- 
perature has been observed, and the idea of using shot 
noise measurements to learn more about potential de- 
phasing/decoherence mechanisms has been introduced. 
Recent experimental results present a puzzling picture 
(e.g. oscillations in the visibility^^), that has not been ex- 
plained so far. Under the assumption that at least part of 
the loss in visibility is due to decoherence processes, the 
observed decrease in visibility with increasing bias volt- 
age is a good indication that Pauli blocking effects are 
important, as this effect is due to lifting the restrictions 
of Pauli blocking on the scattering of particles. 

On the theoretical side, the loss of interference contrast 
in the current had been studied for the Mach-Zehnder 
setup ^* prior to this experiment. More recently the in- 
fluence of decoherence on shot noise has been analyzediS 
(seei^ for related work in quantum point contact), re- 
vealing important differences between phenomenological 
and microscopic approaches, leading to renewed inves- 
tigations on the so-called dephasing terminal modeU^ 
and calculations of the full counting statistics in the MZ 
setupii. However, all of these works deal with a classical 
noise field acting on the electrons. Thus, experimentally 



(i^t(o,o)7/;„(0,t))= / (dfc) /„fee-»-'=* (11) 

(note h = 1), with a band-cutoff k^. Here and in the 
following, we use the notation (dfc) = dk /{2tt). 

The particles are assumed to have no intrinsic inter- 
action, but are subject to an external free bosonic quan- 
tum field V (linear bath) during their passage through 

the arms L,R: Hmt = J2\=l rI '^■''^>^(^)''^^(-'^) '^^ith 
nx{x) = ipl{x)ipx{x). 

We focus on the current going into output port 3, 
which is related to the density: I{t) = evpn3{t) with 

"-3(0 = '03t'03t J where we take fields if^at = V'a(^j t) at the 
position of the final beamsplitter B (except where noted 
otherwise). In the following we set e = vp = 1, except 
where needed for clarity. We thus have 



= Rb [i^'L^L)+TB ^^jjVi^^+e^^i^rs ^Vj^V-L^+cc. 

(12) 
Therefore, the calculation of the average current has been 
reduced to a calculation of the elements of a density ma- 
trix ( il}y ijjx > describing the coherence properties of the 

fermions right at the second beam splitter (after having 
been subject to the quantum noise field). We have set 
Tb = \tB? and Rb = \- Tb- 



C. Influence on the interference contrast 

In this section, we will remind the reader of our results 
for the influence of the quantum bath on the interference 
term in the current I{(f>). These have already been pre- 
sented in a brief communication^, but we repeat them 
here in order to keep the discussion self-contained. They 
form the basis of the subsequent sections on the current 
noise. 

In order to obtain the interference term in the cur- 
rent, we expand the exponential, Eq. (O, to second or- 
der, insert the formal solution, Eq. ©, and perform 
Wick's averaging over fermion fields, while implementing 
a "Golden Rule approximation" , i.e. keeping only terms 
linear in the time-of-flight r. 

These steps will be explained in more detail below, in 
Section FlV CI for the case of the current noise, so we do 
not display them here. 

Note that accounting for cross-correlations between 
the fluctuations in both arms ("vertex-corrections") is 
straightforward for a geometry with symmetric coupling 
to parallel arms at a distance d (assuming d <C Z). Then, 
in the following results (Eqs. lO, (|nj-(I33, and r,^), 

we have to set (vv\ = (VlVl) - (VlVr) and Z?^ = 

D^j^ — D^j^. These correlators, of fields being defined on 
the one-dimensional interferometer arms, actually have 
to be derived from their threcdimensional versions, e.g. 

lvL{x,t)VR{x',t)) = iv{x,y + d,z,t)V{x',y,z,t')) if 

the arms are parallel to the x-axis and separated in the 
y-direction. 

Without bath, the interference term is given by 



^r^l) ^ = rAt\ / {dk)5fk = rAt\{eVl2Ti), (13) 

where we define 5fk = Jik - f2k and fk = {fik + ./2A;)/2 
(for later). 

The leading correction to the interference term can be 
expressed in terms of a phase-shift and a dephasing rate: 



Hi'l^L 



Note that the "classical" contributions (ipl^ij^x 

A = L, R) are not affected by the noisy environment. 
Here the effective average phase shift induced by cou- 
pling to the bath is energy-dependent, and given by: 



rAt*A / {dk)Sfk[iS^{k) - r^(fc)T] (14) 



(with 



6<p{k) =. t{Ra - Ta) f{dq)(ReD 



R 
9,9 



Ko)Sh- 



(15) 



Essentially, the phase shift is due to the effective coupling 
between the electrons, mediated by the bath (containing 
Hartree and Fock contributions). For that reason, it de- 
pends on the nonequilibrium Fermi distribution (differ- 
ence) 5f. The phase shift vanishes for Ta = 1/2, since 
then there is complete symmetry between both arms. 



The suppression of the interference term is quantified 
by the dephasing rate r^(fc), within the Markoff/Golden 
Rule approximation adopted here. In the case of classical 
Gaussian noise, the suppression can be evaluated exactly 
("to all orders" in the system-bath interaction). It is 
equal to exp(— (<y3^) /2), where ip is the phase difference 
between the two arms of the interferometer, fiuctuating 
due to the action of the noisy potential. For the case of 
a single particle coupled to a quantum bath, the same 
suppression factor would be given in general by the over- 
lap of bath states that have evolved under the infiuence 
of the particle traveling along the left or the right armi^. 
Up to now, we have not been able to find an equally 
simple interpretation for the many-particle case. 

The total dephasing rate is T^{k) = r^(fc) -I- T^{k). 
For equal coupling to both arms, this can be written as: 

r^(e) = / — D0S,(u;)[2n(a;) + l-(/(e-a;)-/(e+w))] 

Jo VF 

(16) 
The rate (at energy e = e{k) = vpik — kp) + £f) is an in- 
tegral over all possible energy transfers w from and to the 
bath (which have been combined, so w > here). They 
are weighted by the bath spectral "density of states" 
DOSg(w) = —lmD^{uj)/Tr, where q = uj/vp for ballis- 
tic motion (in this definition, DOS has the dimensions 
uj/q). 

The first term in brackets, 2n{u!) + 1, describes the 
strength of thermal and quantum fiuctuations (with 
n{Lu) = {e^/"^ — 1)^^ the Bose-Einstein distribution). It 
stems from the V(o) in the quantum phase. By itself, 
this would give rise to an energy-independent rate and 
a visibility independent of bias voltage, in contradiction 
to experimental results. In fact, such a procedure (drop- 
ping the back-action terms) would describe a different 
physical situation: that of a single particle coupled to a 
quantum bath (in absence of the Fermi sea) . 

Thus, the second term is crucially important. The 
"back-action" oc D^fi introduces the nonequilibrium 
Fermi functions (/l = Ra/i + 1^/2, Ir = Ta/i + i?A/2, 
and their average, / = {h + /i?)/2^ = (/i + /2)/2) 
which capture the physics of Pauli blocking: Large energy 
transfers vp\q\ ^ eV,T are forbidden for states k within 
the transport voltage window. This can be seen in Fig. |31 
which displays the energy-dependence of the dephasing 
rate, as a function of voltage and temperature. For the 
simplest example of an optical phonon mode (where only 
an energy transfer ujq is allowed), we find two dips in the 
dephasing rate at large voltages. These occur around the 
edges of the non-equilibrium Fermi distribution f, i.e. at 
the edges of the voltage window, and their width is 2ujo. 
When the voltage is reduced, these two dips merge and 
the rate goes down to zero. Thus, when averaging this 
rate over the voltage window (in which electrons con- 
tribute to the current), the average rate becomes zero for 
V,T ^ 0. As a result, the interference contrast (visi- 
bility) becomes perfect (see also the energy-averaged de- 
phasing rate depicted in^). In contrast, at higher tem- 
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Figure 3: Energy-resolved dephasing rate for an optical 
phonon mode (at wo), as a function of transport voltage ap- 
plied to the Mach-Zehnder, for two different temperatures: 
(a) r = 0.1a;o, (b) T = O.Si^o- Here To = r^(e -» oo,T = 0). 



perature, two effects increase the dephasing rate: First, 
thermal smearing of the Fermi distributions reduces the 
restrictions of PauU blocking, and second, thermal fluc- 
tuations in the bath lead to processes of induced emission 
and absorption. 

Note that the strong energy-dependence of the dephas- 
ing rate in the many-fermion case is markedly different 
from the single-particle situation, and thus the depen- 
dence on the bath spectrum is completely different as 
well. In the single-particle case, it is enough to know the 
variance \f/ of the fluctuating phase difference, in order 
to calculate the loss of visibility. In the many-particle 
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Figure 4: Contributions to the decoherence rate in a Keldysh- 
diagrammatic treatment. Left: Diagram involving both ther- 
mal and quantum fluctuations of the bath, but no Fermi dis- 
tributions. This diagram is the same in a single-particle sit- 
uation. Right: Diagram corresponding to the "back-action" 
term discussed in the equations of motion approach. It in- 
volves the fermionic Keldysh Green's function (that contains 
the Fermi distribution, G'^ {e — uj) ocl — 2/(e — uj)) and the 
bath's retarded propagator (describing the response). 



case, we have to keep track of the full bath spectrum 

{vv) . 

As we have only evaluated the corrections to lowest or- 
der, we should be able to make contact to Fermi's Golden 
Rule, describing the scattering of electrons inside the in- 
terferometer arms, by emission or absorption of phonons 
(bath quanta). Indeed, it turns out that the dephasing 
rate is related to Golden Rule scattering rates. How- 
ever, we emphasize that it is not given solely by the rate 
for scattering of particles, as one might naively assume. 
Rather, hole-scattering processes provide an equally im- 
portant contribution to the dephasing rate, which is thus 
the sum of particle- and hole-scattering rates. In our 
case, we find: 






^ p 



L/R 



^L/R 



)/2, 



(17) 



with r 
rL/R. 



L/R 



(k) - Jidq) (vv) (1 - h/R.k-q) and 
\ ' q-q 

„ '(k) = jidq) (yv) fL/R.k+q- This is because 

\ '9,9 

both processes destroy the superposition of many-particle 
states that is created when a particle passes through the 
first beam splitter, entering the left or the right arm. A 
more detailed qualitative discussion may be found iniS, 
for the case of weak localization, and in the next subsec- 
tion. 

For linear transport, i.e. a the limit of infinitesi- 
mal bias voltage F ^ 0, we have fhk-q — fhk+q -^ 
— tanh(/3(A; — g)/2) under the integral. Then we recover 
the result well known in the theory of weak localizations^, 
where ballistic motion in our case {uj = vpq) is replaced 
by diffusion. 

Finally, we note that a treatment using Keldysh dia- 
grams would yield (in the absence of vertex corrections) 
a dephasing rate that is equal to the decay rate of the 
retarded (or advanced) propagator, and thus given by 
the sum of the two diagrams shown in Fig. ^ These 
correspond exactly to the first and the second contribu- 
tion discussed above. For the average current, the effort 
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Figure 5: Contribution of particle- and hole-scattering pro- 
cesses to the dephasing rate in a many-fermion interferome- 
ter. 



coherent superposition that requires not only the pres- 
ence of a particle in one arm but also the absence of a par- 
ticle in the other arm. This is why the many-body super- 
position can equally be destroyed by particle- and hole- 
scattering (leading to states with | . . . , 0, . . . ; . . . , 0, . . .) or 
|. . . , 1, . . . ; . . . , 1, . . .), respectively). We emphasize that 
the dephasing rate is independent of the amplitudes Ia 
and TA in this superposition. The reason is basically that 
the dephasing rate describes the decay of the off-diagonal 
element of the density matrix (in the space of these two 
states) , and that the amplitudes only enter as a constant 
prefactor in that element. Thus, the dephasing rate is 
simply given by the sum of particle- and hole-scattering 
rates, as noted above. The factor 1/2 arises because we 
are not asking about the decay of populations (which is 
described by Tp and Th) but essentially the decay of a 
wave amplitude. This is the same factor that arises in the 
relation T2 = 2Ti known for pure dephasing processes in 
the context of Bloch equations. 



involved in both calculations (Keldysh or equations of 
motion) is still about the same (a few lines). However, 
for the shot noise corrections discussed below, we found 
the equations of motion method much more convenient. 



D. Particle- and hole-scattering contributions to 
the dephasing rate 

In this section, we briefly provide a more qualitative 
discussion of the fact that hole-scattering processes lead 
to an equally important contribution to the dephasing 
rate F^ = (Fp -|- F;i)/2. The ratio of Tp and Th depends 
on the energy under consideration, with Fp providing the 
full dephasing rate at high energies, and F^ accounting 
for F^ at low energies (seeS). 

This is a generic feature for decoherence of fermionic 
systems. Even though it is implicit in known diagram- 
matic results^, we are not aware of any simple physical 
discussion (except our own recent treatmenljiS in the case 
of weak localization). From the perspective of a single 
particle, the first beam splitter creates a superposition of 
the form tA\R) +rA\L), with the states R/L denoting a 
wave packet inside the right/left arm. In the presence of 
a sea of other fermions inside the interferometer arms, we 
should write instead a superposition of many-body states 
(see Fig. |5j, schematically: 



tA\ 



,Q, 



.,l,...)+rA|...,l,. 



,Q,...) 



(18) 

We have indicated the occupations |left; right) of single- 
particle states in both arms, with a bar denoting the en- 
ergy level e of interest and the remaining particles (in the 
nonequilibrium distributions) playing the role of specta- 
tors. The interference term (i'l'ipRj is sensitive to the 



IV. CURRENT NOISE IN THE 
MACH-ZEHNDER SETUP 

A. Introduction 

As our method yields directly the modified particle 
fields, it may be used, in principle, to calculate any 
higher-order correlator of those fields. Of particular ex- 
perimental interest is the current noise in the output port 
of the interferometer. This has been (and is) currently 
being studied in the Weizmann MZ setup^^^ 



B. General properties 

The zero-frequency current noise power is defined as 



S 



+ 00 



di((/»/(0) 



(19) 



where the double bracket denotes the irreducible part: 

/\ . For the MZ setup 

considered here, the current noise only has contributions 
up to the second harmonic in the external flux: 



mm)) = (mm 



S = So + Si cos{(j) - S(j)i) + S2 cos(2(0 - 5<P2)). (20) 
The dependence on cf) and Tg, Rb can be made explicit, 

S ~ RbTbCq + RbCor. + TgCoT + (21) 

2Re [e'^(t|jrB)(i?BCifl + TbCit) - e^">'TBRBC2\ . 

with the coefficients following directly from inserting Eq. 
into lfTn|l . see below. Here S'o, 5*1, S'2, (50i and (502 can 
be obtained by comparing Eqs. H21|) and H2U|I . 



The coefficients Co, Cqr, . . . are expressed in terms of 
four-point Green's functions, similar to the expression 
for the average current. These, in turn, contain the fuh 
dependence on interactions, as well as on voltage, tem- 
perature, and Ta- We list them for reference, setting 
Lt = ijJL{t) and Rt = 'ipR{t) for brevity. We will also set 
wf = 1, as before. 



under cj) t—f —(f> and V ^-^ —V separately, but these sym- 
metries may be broken by a bath-induced phase-shift, to 
be discussed below. 



C. Evaluation of current noise to leading order in 
the interaction 



Cor = I dt (^(LlLtLlLo 



Cot = I dt{{ RlRtRlRo 



(22) 
(23) 



Co = I dt{{ LlLtRlRo + RJRtLlLo ) ) + (24) 



dt {{LlRtRlLo)) + ((RlLtLlRo)){25) 



CiR = j dt (^(^R\LtLlLo)) + (^(^L\LtRlLo)) {2&) 
CiT - f dt (iRlRtRlLo)) + UR\LtRlRo))m 



C2 - I dt ^{rIurILo 



(28) 



Cqir/t) are real-valued, the other coefficients may be- 
come complex. 

In the absence of a quantum bath, these coefficients 
have the following values: 



C, 



(0) 
OR/T 



(dk) 



h{l~h)~-^{RA-TAf5fl 



(29) 



^0 



cr = i[dk){fLk[i-fR.k) + iRk{i-fLk))}m 



^2RaTa [(dk) 5fl (31) 

Cf^)/j, = ±rAt*A{TA-RA) j{dk)5fl (32) 

Cf = RATAJ{dk)5fl. (33) 



Those expressions yield the result given by the well- 
known scattering theory of shot noise of non-interacting 
fermion s'^^''^^-^^'^'^ : 



In order to evaluate the correlators lf^ - (PS|) to leading 
order in the interaction, we expand the general solution 
of the equations of motion for the electron operators. Let 
Lj denote the unperturbed electron field, and g a formal 
expansion parameter (to be set to 1 in the end). Then 
we have, for the electron field at the end of the left arm, 
just before the final beamsplitter: 



U 



1-ig dtiVL{ti,t)- 



(35) 
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ti 



dti / dt2VLiti,t)VL{t2,t) 

Jo 

tg' I dti I dt2 I dx2 X 
Jo 



-(o)t 



(0) 



^LL(^l:*;2;2,t2)^^+J2_^^/^j,i^+t^_^^/„^ 



(0) 



We have expressed the arguments of the potentials V 
and response kernel D^ in terms of the time ti = . . . t 
elapsed since entry into the left interferometer arm, with 
the electron moving from x = to 2: = ^ = vpT during 
the corresponding time- interval [t — T^t]. We have set 



VL{tl,t) = VL{vFtut^T + ti) (36) 

L)fi(il,i;2:2,i2) = D^{vFtl-X2,t-T + ti-t2lV) 

assuming a stationary environment that is translation- 
ally invariant. The expressions for Rt are completely 
analogous. In writing down Eq. (|35|l . we have omit- 
ted the cross-term Df^j^, assuming that the wavelength 
of relevant fiuctuations is considerably shorter than the 
distance between the arms of the interferometer (such a 
term can be added easily, see the remark above, in Sec- 
tion I^^P). This also implies (VlVr) ^ 0. In terms of 
the bath spectra, we have (both for L and R): 



'(0) 



(dfc)(/2fe + SfkT){l - {hk + SfkT)) , (34) 



where r(0) = TaTr + RaRb + 2t\rAt*BrB cos(0) is the 
transmission probability from 1 to 3. 

For our model, the full shot noise power S may be 
shown to be invariant under each of the following trans- 
formations, if the bath couples equally to both arms 
of the interferometer: (i) ^a ^^ tAt 4> ^~^ ~4' (ii) 
V I— > — y, 4> >— > —(f) (iii) ts '^ tb- As a consequence, 
CiT = —CiR. Note that the free result H34|l is invariant 



VL{t'l,t')VL{tl,t)) = 

(dq) j {duj)e'^i^^<i--')(^'i~*^)--'(t'~^)Uvv\ (38) 

bf^{ti,t-X2M) = 

(dq) /"(da;)e*[9(^i=-*i-^2)-'^(*i-*2+t-r)]^fl (39^ 

We now evaluate the leading order (g^) correction to 
the noise power l|21|) . by inserting the expressions for Lt 
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and Rt into the coefBcients Co,CiT,C'ifl, and C2 (Eqs. 
(I13-(|2Hl)- Bare electron operators are contracted using 
Wick's theorem, and the resulting averages can be per- 
formed by expressing L[ ' , R\ ' via ^1,2 (Eqs. Q,®) 



and employing Eq. lfTT|l . 
-" — -"7 \ 



resentations (VV) and D^ 



After inserting the Fourier rep- 
all temporal and spatial 

quj 

integrations have to be carried out. In doing so, we will 
use a Golden Rule (Markoff ) approximation, i.e. we keep 
only the leading order in r. 



dti 



du^iMti^t^) 



X + iO 



(40) 



(and so on), assuming the correlation time of the en- 
vironment to be much shorter than the time-of-flight 
T. Although it is in principle straightforward to go be- 
yond this approximation (evaluating all these integrals 
exactly), the result gets very unwieldy, and other effects 
(such as the curvature of the interferometer paths) should 
be taken into account as well on that refined level of de- 
scription. Thus, we are neglecting the fact that energy- 
and momentum-conservation will only be fulfilled up to a 
Heisenberg uncertainty t~^ and l~^, respectively. Within 
this approximation, we have been allowed to extend the 
X2-integral in Eq. (|35|l over all of space, even though the 
interaction is assumed to be restricted to the interferom- 
eter arm (it will be restricted automatically by the short 
range of D^ and the fact that ti G [0, t]). 



D. Current noise corrections due to the quantum 
bath 



SClR 



{dk){dq)ReD^,j x 



Im 

[-SfkSfk+q{6fk+q + 2dfk){TA - Ra)V2+ 
2RATA6fk5fl+^ + Sfk+Jki3 - 2/fc) 

-^fkfk+q + 2fkfk+q{6fk - Sfk+q)] 

+D"'fi{GV/2TT) [{dk) X 



Sfk(.^{Tl + R\) - 5RaTa) - 2A(1 - h) 



2tRaTa ^ ' ' \ /o,o 

2 f{dk)idq) ImD^JkSfk+qiSfk + Sfk+q) 



,(43) 



idk)idq){VV) 5fk{5fk + 5fk+q) 
I q-q 



Im<5C2 - -ARaTa / {dk)5fl5^k 



(44) 



(45) 



To obtain physical insights, it is best to translate the 
coefficients that have been obtained above into correc- 
tions to the different harmonics 5*0 , 5*1 , S2 of the noise 
pattern S{<j)) and the phase shifts 5(j)i and 5(j)2 (compare 
Eqs. (pn|l and (|^ ). Then we find, for the lowest order 
corrections: 



After a straightforward but lengthy calculation, we ar- 
rive at the leading-order corrections to the coefficients 
Co, Cm, Ci in the noise power S . Here we list the explicit 
analytical results for the shot noise correction (cf. Eq. 
I21|l . valid for arbitrary bath spectra (note 5Cqb.it — 
and 5CiT = —SCm): 



SCo 



idk)idq)lmD^;x 



^tRaTa 

[SfkSfk+qifk+q - fk) + 

iflk + f2k)fk+q — iflk+q + l2k+q)Ik\ + 

(dk){dq) (vv) X 

[iflk+q - /lfe)(l - /ife) + i.f2k+q - ./2fc)(l - f2k)] 



\ /o,o 



(41) 



Re 



(5C- 



\R 



TrAt*A{RA-TA) 



{dk){dq)lu\Di 



[SfkSfk + q{fk+q + 3/fe - 2) + Sf, + Jk - 6flfk+q] + 



{dq){VV) 

I q,q 



{dk)5fl 



(42) 



^^0 
5Si 

5S^ 
S2 



RbTbSCq 
ilmSCiR 

w 



ReSC2 

~^ 
iReSCiB, 

MO) 

lmSC2 

~~C^ 



(46) 
(47) 

(48) 

(49) 

(50) 



(Note, when comparing with Eqs. (|42|l and (|43|l . that we 

took into account rAt*A, and thus also C-J^, being purely 
imaginary) 



E. Discussion of current noise in the 
Mach-Zehnder coupled to a quantum bath 

The results of evaluating Eqs. (|46|I - H50|I are shown in 
Figs. © and Q for the illustrative example of a damped 
optical phonon mode, D^^ = a[{uj — wq + if])^^ — {uj + 
1^0 + ''■v)~^]j with ri/uo — 0.1. 



<5S2/(c<tS2) 




(a) -fw^ct>i/(ar) (b) S<p2/((xr) (c) 



Figure 6: (a) Correction to the flux-averaged current noise 
power So for a damped optical phonon mode (strength a, 
time-of- flight r), normalized with respect to unperturbed 
value, (b) and (c): The corrections to the first and second 
harmonics Si and S2 are negative, revealing the loss of inter- 
ference contrast in S((/f>). Energies are plotted in units of the 
phonon mode {ujo = 1), and the MZ setup has been chosen 
asymmetric {Ta = 0.3, Tb = 0.4). 




L T 



eV ' ^-1 




Figure 7: Corrections to the phase shifts 5(j>i (a) and 5(j)2 (b) of 
the first and second harmonic in the shot noise pattern 5(0), 
in an asymmetric MZ setup (parameters as in Fig. (SJ. These 
phase shifts will, in general, be different from the effective 
phase shift Sip in the current pattern, leading to the situation 
schematically depicted in (c). 



As expected, the (/)-dependence of the shot noise H21|l is 
suppressed, i.e. not only the visibility (interference con- 
trast) of the current pattern /(^) but also that of the 
shot noise pattern S{(j)) is reduced by the bath: see Fig. 
|Sl(b) and (c). We emphasize that this reduction becomes 
noticeable only once the voltage or the temperature be- 
come comparable to the frequency of the phonon mode. 
Only then the particle can lose its coherence by leaving 
a trace in the bath (that acts as a kind of "which-way 
detector"). This is the same behaviour found for the 
visibility of the current, and it is satisfying that this sim- 
ple qualitative physical idea also holds for dccohcrcncc in 
shot noise. Note, however, that we have not found a way 
to express the comparatively complicated formulas for 
5Si and SS2 in terms of the simple dephasing rate which 
we derived above, Eq. ^] It is interesting to note that 
the decrease of the second harmonic S2 proceeds faster 
than that of the first harmonic. Si. This is qualitatively 
consistent with the observations made by Chung et al. 
for a MZ setup using the phenomenological dephasing 
terminal model^^. 

There is no Nyquist noise correction, as seen in Fig. 
El (a), at V — 0. This can be understood easily, since 
the (unperturbed) Nyquist noise 5(0) (1^ = 0) does not 
depend on </> and thus should not be sensitive to a noisy 
environment that changes the phase (j). 

The limit of classical noise (treated to all orders in 
Refsi^) is recovered by setting D^ = and using the 

symmetrized correlator (K;iK:i) — (<V,V>)/2 every- 
where in the shot noise correction derived here, with the 
exception of Eq. H41|l . which has to be replaced by: 



scs' 



{dk)idq){V,,Vcl)^^x 



[(/ife+g - /Lfe)(l - fBk) + [fRk+q - .fRk){i - fLk)] + 



'^RaTa 



[dk)5fk 



(^clKl)o^O 



(51) 



This contribution contains a finite ^-independent 



Nyquist noise correction (cf.^*'), in contrast to our result 
for the quantum bath. This may be understood as be- 
ing due to heating of the MZ electrons by a bath which 
is nominally at infinite temperature (according to the 
fiuctuation-dissipation theorem FDT, applied to the case 
i^^ = 0). 

We emphasize that it is impossible to recover the full 
quantum noise result by inserting some suitably modified 
classical noise correlator {Vc\Vc\). This is in contrast to 
the dephasing rate, where such a procedure (with {Vc\Vc\) 
containing Fermi functions for Pauli blocking, see Refs.^^ 
or also2&) can be made to work. In particular, having 
only classical noise cannot yield the important phase shift 
terms. In contrast, the conductance fluctuations are cor- 
rectly captured even by the classical approach. 

At large V (larger than the bath spectrum cutoff), 
there is a contribution oc V'^ in 5Sq and 882, due to time- 
dependent conductance fluctuations {51 oc SG{(f){t)) ■ V), 
corresponding to the leading order of "Sci" in Refs^ (see 

(vv\ terms in Eqs. (gTJ and ll^ ). 

As mentioned in^ the main surprising feature con- 
nected to the shot noise correction is the behaviour of 
the phase-shifts (50i.2- Naively, one might expect the 
effective phase shift to be one and the same for all quan- 
tities depending on the Aharonov-Bohm phase, whether 
it be the current /((/)) or the shot noise S{(f>). However, 
the phase-shift S(j)2 in the e^"'^ term is twice as large as 
expected from the phase-shift in /(^) and, moreover, the 
e"^ phase-shift S(pi does not vanish even if Ta = 1/2 
(but Tb ^ 1/2). As a consequence, and in contrast to 
the current /(</>), even for completely symmetric inter- 
ferometer arms (same density, same Fermi distributions, 
same coupling to the bath), there remains a <-> — 0- 
asymmetry in bS. The explanation* rests on the fact that 
the phase shift is sensitive to the density difference be- 
tween the arms (as discussed above). As a consequence, 
density fluctuations in both arms also lead to fluctuations 
of this phase shift. While the average current only feels 
the average phase shift 8(p^ the current noise is affected 
by those fluctuations. The extra terms in (55, which are 
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responsible for the deviation from the behaviour of the 
average current, come about because the fluctuations of 
the phase shift are correlated with the output current, 

{dip{t) - Sif){i{0) - I)) ^ 0. This is a straightforward 



consequence of the fact that the output current / itself 
is correlated with the currents/densities traveling inside 
the interferometer arms. This also explains the fact that 
Ta ~ 1/2 is not enough to obtain a (/)-symmetric shot 

noise (since the correlator / 8(pl ) depends on Tb as well) . 

We emphasize that a fluctuating effective phase shift 
depending on the density fluctuations inside the arms 
will quite likely be present in any model of interacting 
fermions moving inside an interferometer (either with in- 
trinsic interactions, i.e. as a Luttinger liquid^, or with 
interactions mediated by a bath, like in the present work). 
Thus, the consequences (different phase shifts in current 
and current noise, and different phases of the two har- 
monics in the current noise) will hold more generally than 
our specific model. It is thus important to carry out ex- 
periments that test for those phase shifts in asymmetric 
interacting interferometers. 



ing many particles coupled to a quantum bath. It takes 
into account the simplifications provided by the chiral 
motion at approximately constant velocity, and is thus 
more efhcient than more general approaches. In particu- 
lar, it is able to keep, in a straightforward and physically 
transparent manner, many-body effects, such as Pauli 
blocking (described as a consequence of the backaction 
of the bath onto the system) or the influence of hole- 
scattering processes in the case of fermions. We have 
applied this method to the fermionic Mach-Zehnder in- 
terferometer, presenting full analytical results for the in- 
fluence of the quantum bath on the current noise. As we 
have discussed, the main effects are a reduction of the 
interference contrast in the shot noise pattern S{(j)) and 
a peculiar behaviour of the effective phase shifts in the 
two harmonics of S{(j)), for asymmetric setups. 

We are anticipating future applications such as the 
treatment of higher-order effects of the bath or decohcr- 
ence in bosonic (atom-chip) interferometers. 
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